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Abstract 

We study non conformal negative modes of Schwarzschild black holes in the five di- 
mensional Einstein-Gauss-Bonnet gravity theory with a zero or negative cosmological 
constant, in the context of the semiclassical quantum gravity formulated by the path 
integral. We find an unique negative mode when the black hole has a negative heat 
capacity which is the same in Einstein theory. On the other hand, we still find negative 
modes inside a parameter region where a small black hole has a positive heat capacity. 
The number is one/two for the case of zero/negative cosmological constant respectively. 
In the rest of the parameter region where the heat capacity is positive, we find no nega- 
tive modes. We discuss the possible physical understanding of having one or two negative 
modes. 



1 Introduction 



Black holes are very important objects in classical and quantum gravity. One of the most 
important and profound properties in black hole physics is its thermodynamics. The thermo- 
dynamic quantities can be derived from the path integral formulation of quantum gravity, and 
thus have connections to quantum nature of gravity. 

In the case of four dimensional Euclidean Einstein quantum gravity, the Schwarzschild 
black hole is an instanton to induce a decay of hot flat space by nucleation of black holes. 
This instability is suggested by a unique off-shell negative mode around the Schwarzschild 
instanton p.J. The existence of negative mode is believed to have a connection with the 
local thermodynamic stability of Schwarzschild black hole. For example, the Anti de-Sitter 
(AdS) Schwarzschild black hole has an unique negative mode when the black hole has a 
negative heat capacity, and the negative mode disappears exactly when the black hole becomes 
local thermodynamic stable [2]. Reall pj discussed in detail the connection between the 
local thermodynamic stability and the stability against small perturbations around a black 
hole solution (existence of negative modes) and showed they agree under some reasonable 
assumptions. 

The negative mode has been also studied in the connection with Gregory-Laflamme insta- 
bility [Ij, since the negative eigen value is related with a Kaluza-Klein mass square of graviton 
along the directions to which a black string extends [3]. Gubser and Mitra [5] conjectured that 
there is a negative mode iff the black hole has a thermodynamic instability (the correlated 
stability conjecture). This conjecture is known to be satisfied for many different black branes 
including the black branes discussed by Gregory and Laflamme [1] , black p-brane solutions in 
string theory [6] and also black strings in Anti de-Sitter space [21 El [7]. Also more complicated 
systems such as the D0-D2 bound state [8] and the non-extremal smeared black branes [9] are 
discussed and the conjecture is satisfied. There are, on the other hand, counter examples [TO] , 
and those are related with the non uniqueness of black string under fixed conserved charges. 
The correlated stability conjecture is checked to be true for all the known cases when a black 
string is uniquely determined by the conserved charges. 

Since the black hole thermodynamics has a deep connection with quantum mechanics of 
gravity, we have a natural question whether the properties related with the black hole ther- 
modynamics are modified, if we add higher dimensional corrections to Einstein Hilbert action 
which are expected to appear in an effective theory of quantum gravity, for example, string 
theory. The Gauss-Bonnet term is one of such terms and is the lowest dimensional term among 
them. The black hole solutions are known [Tl] and the thermodynamics is also calculated [T2] 
in Einstein- Gauss-Bonnet theory. We naively expect that the local thermodynamic stability 
and the stability against small perturbations still agrees in Einstein- Gauss-Bonnet theory. In 
this paper we study the negative mode of Schwarzschild black hole in Gauss-Bonnet theory 
with and without a negative cosmological constant. Since the Gauss-Bonnet term becomes 
topological in four dimensions, we study the five dimensions. 
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We find there is an unique negative mode when the black hole has a negative heat capacity. 
As the black hole becomes larger, the negative eigenvalue goes up and becomes zero exactly 
when the sign of heat capacity changes from negative to positive. This behavior is exactly 
the same as that in Einstein gravity and consistent with the fact that the effect from Gauss- 
Bonnet term is negligible when a black hole has a large size and the spacetime is weakly 
curved. However as the five dimensional black hole becomes smaller and starts having a 
positive heat capacity due to the Gauss-Bonnet term, we find a new behavior. We still 
find negative modes and the numbers of negative modes are one or two for the case of zero 
or negative cosmological constant respectively. We also find one or two positive modes for 
the case of negative cosmological constant. As the black hole becomes further smaller, the 
negative modes disappear. These results lead us to reconsider the connection between the 
thermodynamic stability and the existence of negative modes. 

In the next section, we review the thermodynamics of the Schwarzschild black hole in 
Gauss-Bonnet theory. Then in section 3, we carry on the S-mode perturbations around the 
black hole and study the non conformal negative modes. We only study a simple case, i.e. 
the five dimensional Gauss-Bonnet theory. The detail procedure of numerical analysis is 
summarized in Appendix A. We summarize our numerical results as figures in FigHJIS] (and 
also Figl3l8]in Appendix B). The real numerical data is listed in Appendix B. In section 4 we 
interpret our numerical results and conclude in section 5. 



2 Black hole thermodynamics 

In this section we summarize the setup and thermodynamic quantities of black hole in Gauss- 
Bonnet theory. 

The path integral formalism of Euclidean quantum gravity with a finite temperature T is 
schematically written 

Z = Jvgexp[-I{g)], F = -T\nZ, (2.1) 
IGnGn Jm 

J-'GB — K — ^Kabri + Uabcd-ti , (,^.Oj 

where F is the free energy, Gn is the Newton constant, A is the cosmological constant, c is a 
coupling constant of mass dimension —2 in front of Gauss-Bonnet term Lgb, h is the induced 
metric on the boundary, 7^ is a surface term and n is the spacetime dimensions. We also 
need gauge fixing terms. In this section we consider A < 0. The manifold is asymptotically 

X and the periodicity for along the Euclidean time r is given t = t + f3, f3 = 1/T. 

The partition function Z is semiclassically approximated by the sum of classical paths 
(saddle point approximation). We are then interested in a black hole solution and the small 
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perturbations about the solution. The black hole solution in the n dimensions is known [TT] 



ds^ = Qnhdx^'dx^ = f(r)dT^ + + r^dQl_2, 
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where is a parameter and is related with the black hole mass M = yU^„_2/(87rG„,) where 
An-2 is the area of unit {n — 2) sphere. There are two solutions e = ± and the solution with 
e = — 1 becomes the Schwarzschild black hole in Einstein gravity in the limit of c — > 0. Thus 
we only discuss the case of e = — 1 in this paper. The temperature T is determined to avoid 
a conical singularity at the horizon and then given by the surface gravitational force at the 
horizon r = rt, (/(r^) = 0), 



^ ^ f'jrh) ^ {n-2){n- 3) {2r^ + c(n - 4)(n - 5)} - 4r^A 
47r Sirrhin - 2) {rl + c{n - 3){n - A)} 

for > — c(n — 3)(n — 4). Using the horizon radius r^, the black hole mass is written 



(2.6) 
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The entropy which satisfies the first law in thermodynamic mechanics is computed 
I{g) = ^M{g)-S{g), 



S 



A 



4G. 



[rr^ + n{n - l)crr^] 



(2.7) 
from 

(2.8) 



where the entropy is no longer simply proportional to the area of the horizon. The local 
stability as a thermodynamic system is determined from the sign of heat capacity Cy, 



dM 
'df' 



(2.9) 



If we denote a one parameter family of geometries which contains the black hole solution as 
X, the second derivative of the action along the path x is 



dx'^ 
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(2.10) 
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Thus if the specific heat is negative, the perturbations around the black hole where {dI/dx)T = 
must have at least one negative mode [3] . This is a basic observation that the local thermo- 
djTiamic stability is equivalent with the stability against the perturbations about the black 
hole solution in the semi classical quantum gravity. 



3 



Black hole in five dimensions 



The heat capacity depends on various parameters. Hereafter we study n = 5, since the five 
dimensions is the lowest dimension where the Gauss-Bonnet term non trivially contributes 
to the local dynamics of gravity. In five dimensions {n = 5), the mass is a monotonically 
increasing function in terms of Vh and then the sign of heat capacity is equal to the sign of 
dT/drh- Since the temperature is 

T = ^~^h^ (2 11) 

67r(r,+ (2cK))' ^ ' ' 

for rl > — c/2, we divide three regions (1) c < 0, (2) < c < and (3) c > depending 
on the behavior of heat capacity in terms of the black hole radius Vh- 

(1) c < : The temperature behaves similar to the case of AdS Schwarzschild black hole 
in Einstein gravity, but goes to infinity at Vh = v— 2c before Vh = 0. See FigUl The 
temperature takes the lowest value at Vh = rc, 



^ 3 + 6cA + ^3(1 + 6cA)(3 + 2cA) . (2.12) 



-2A 

Thus Cy < for < and Cy > for > Tc- 

(2) < c < : The temperature is largely modified at a small black hole region and 
there is a region where there are three different black holes at a same temperature. See 
Fig 131 The temperature starts from zero and increases up to the black hole radius r^, 
and starts decreasing until rh = r^. After that, the temperature increases again as the 
black hole becomes larger. Tc is same as (12.121) and is 



^ 3 + 6cA - V3(l + 6cA)(3 + 2cA) , (2.13) 



-2A 



As c ^ ^1^, Ta approaches to i.e. — >• r^. Tc does not change much by changing c. 
Cy < for Ta < Th < T^, and Cy > for Vh < and Vh > r^. 

(3) < c : The region where the black hole has a negative heat capacity disappears and 
the black hole is always stable as a thermodynamic system. See FiglH Cy > for all Vh- 

The behavior of temperature of black hole in five dimensions is special and that in more than 
five dimensions is almost the same as that in Einstein gravity. 



3 S-mode perturbations and negative modes 

We have described the thermodynamic properties of black hole in Gauss-Bonnet theory. We 
now study small perturbations around the black hole solution in five dimensions. The action 
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around the soliton up to the quadratic orders in the perturbations Sgab, {gab — Qab + ^gab), 
becomes 

I = Io(g) + Hg, Sg), h{-g, 5g) = J d^x dg^'AYM^g/, (3.1) 

where the operator A(^) is called the Euclidean Lichnerowicz operator in Einstein theory. 
Thus the one loop quantum collections around the solution are proportional to a/ det A{g) 
and if we have one or more negative eigenvalues in this operator A(^), the black hole solution 
is not a local minimum, but a saddle point. 

We in this paper study the time independent s-mode perturbations around the black 
hole solution. This is because this mode is known to induce the unique negative mode of 
Schwarzschild black hole in Einstein gravity and we expect that negative modes appear only 
in the s-mode perturbations. The time independent s-mode perturbations are given 

SgJ> =dmg{Huir), iJ„,(r), K{r), K{r), K{r)), (3.2) 

where the Lorentz indices are made up and down by the background metric cjab- In the case of 
Einstein gravity, these fluctuations have two independent modes and the one mode is a trace 
mode and the other mode is a traceless mode. The trace mode is related with a conformal 
rescaling and has a ghost like kinetic term. The traceless mode gives the unique negative mode. 
In Gauss-Bonnet theory, there are also two modes and we study the mode which becomes a 
traceless mode in the hmit to Einstein gravity. We name this mode as the 'traceless' mode 
even though it is no longer traceless and study the negative eigenmodes in this mode. Similarly 
we name the other mode which becomes a trace mode in the limit to Einstein gravity as the 
'trace' mode. 

In order to compute the eigenvalues of the operator, we have to solve the equation of 
motion for the eigenmodes 

^VMSg/ = ^Sg%, (3.3) 

with appropriate boundary conditions at the infinity and the horizon. We note that a eigen- 
mode with A 7^ is an off-shell mode and a mode with A = is an on-shell classical per- 
turbation. This eigen equation is complicated and we do not separate the two modes in the 
equation level, but we put the 'trace' mode zero by giving the appropriate boundary condition 
in the numerical analysis. We give the detail procedure of numerical analysis in Appendix A. 

Asymptotic solutions 

We first study the asymptotic behaviors of s-mode perturbations at the horizon and the 
infinity in r. The asymptotic solutions of the 'traceless' mode at the horizon are (e = r — Vh) 

Hu{r) oc 1 + 0(e), //,,(r) oc -1 + 0(e), X«(r) oc ^ + 0{e). (3.4) 

i/«(r) oc - + 0(lne), i/,,(r) oc -- + C»(lne), i^«(r) oc ^ + C)(lne). (3.5) 
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At the leading order, the asymptotic solutions are the same as those in Einstein theory and 
they are traceless. The first solution is a normalizable mode which becomes clear if we use 
Kruskal coordinate. 

Similarly the asymptotic solutions for A 7^ at the infinity are 

Hu{r) oc r^++^ + C(r"+), i^„(r) oc r"+ + 0{r''+~^), Ku{r) oc ^r"++2 + 0(r'^+). (3.6) 
Hu{r) oc r<^-+2 ^ Oir"-), H„{r) oc r"" + 0{r^-^), Ku{r) oc K^-^^ + 0(r'^-). (3.7) 



Here 



a± 



4±2Jl 



A^^ 

c 



2c . 
1 + — A - 1 
3 



(3.8) 



where A is an effective cosmological constant computed from the AdS black hole metric in the 
infinity, i.e. 



fir) 



! + -< 1 

2c 



-I 



^ + A 



(3.9) 



The asymptotic solutions are traceless at the leading order and the second solution (13.71) is a 
normalizable mode. 

Therefore a small perturbation must have the asymptotic solutions (13. 4p and (13.71) at the 
horizon and infinity. In order to obtain the eigenvalues and eigenfunctions, we rely on a 
numerical analysis. In the numerical analysis, we start from a point near the horizon with the 
normalizable mode (13.41) . and solve the equations of motion toward the infinity. We change the 
value of A and find the eigenmode by the shooting method. See more detail in Appendix A. 



Negative modes 

We now show our numerical results of eigenvalue A. Here we only describe the results, and 
in the next section we will discuss how to understand the results. 

As discussed before, the heat capacity depends on the parameter c in front of Gauss-Bonnet 
term and we divide three regions depending on how the heat capacity behaves: (1) c < 0, (2) 
< c < 3^ and (3) c > ^1^. Since the equations of motion are invariant under the following 
rescaling 

r ^ ar, M ^ a'^M, c o?c, A a"^A, A «"^A, (3.10) 

with a > 0, we can always take A = — 1 in the numerical analysis. The numerical data used 
in Fig{l]l5]are given in Appendix B. 

We show our numerical result in FigHlwith (c. A) = (—0.1, —1) in region (1). The behavior 
for a general value c in region (1) is essentially the same. We can see that the behavior is 
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Figure 1: The eigenvalue A in terms of 
for c = —0.1 with A = —1. The value on 
the vertical line is for A (not for T). Vc = 
1.94. A -oo, (T oo), at r/, = 



Figure 2: The eigenvalue A in terms of Vh 
for c = (Einstein theory) with A = — 1. 



similar to the one in Einstein gravity [2l [7] (see Figj2]) except that the eigenvalue goes to 
negative infinity at = y— 2c. We find an unique negative mode when the heat capacity 
is negative. As the horizon size becomes larger, the value becomes larger and goes to zero 
exactly when the heat capacity changes its sign. The eigenvalue is always positive when the 
heat capacity is positive. 

We show our numerical result in FigjS] with (c, A) = (0.01,-1) in region (2). The left 
figure of Figl3] is for a small black hole where the small black hole has a positive heat capacity 
and thus stable as a thermodynamic system. A small black hole with the horizon size € 
[0.024, r^j = 0.143] has two larger black holes with a same temperature. The behavior of 
eigenvalue A is similar to that in region (1) when a black hole has a large size. This is simply 
because a large black hole is weakly curved and the effect of Gauss-Bonnet term is negligible 
in this region. We find no negative mode and one positive mode when a large black hole has 
a positive heat capacity and the eigenvalue decreases as a black hole becomes smaller and 
goes to zero exactly when the heat capacity changes its sign {rh = 1.71). Then the eigenvalue 
becomes negative when a black hole has a negative heat capacity. This behavior is same to the 
previous case. However we find a new behavior in a small black hole case. As the black hole 
becomes more smaller, the heat capacity becomes positive again due to Gauss-Bonnet term 
(^h = 0.143). Although the negative eigenvalue quickly goes to zero, we still find negative 
modes until Vh = 0.094. Moreover we find two negative modes for G [0.094,0.143] and the 
small negative eigenvalue goes to exactly zero at = 0.143. We also find one or two positive 
modes in r/j G [0,0.3]. 
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Figure 3: The eigenvalue A for (c, A) = (0.01, —1) in terms of rh G [0, 0.4] (left) and Th G [0, 6.5] 
(right). Ta = 0.143. Vc = 1.71. A small black hole rh G [0.024, rj has two larger black holes 
with same T. 

We show our numerical results for (c. A) = (1,-1) in FigSlin region (3). In this region the 
heat capacity is always positive and we find no negative mode and only find a positive mode 
whose eigenvalue is essentially constant and does not depend on the black hole size. 

In the final, we show the numerical results for the case of the zero cosmological constant in 
FiglSl We chose c = 0.01. Notice that the equations are invariant under the rescaling (13.101) 
with A = and we can further fix one parameter, say c, without losing generality. The heat 
capacity is positive for small black holes {r^ < V^) on the contrary to the case of Einstein 
gravity where a Schwarzschild black hole always has a negative heat capacity. We find that 
a negative mode still exists even a small black hole has a positive heat capacity as long as 
rh > 0.0612 for c = 0.01. 

We now explain how the behavior of eigenmodes in region (2) continuously reaches to the 
behavior in region (1), (3) or the case for A = 0, by changing c or A. The region (2) approaches 
to the region (1) by decreasing c with A fixed. Since = \/2c+0{c^^'^) for a small c, the small 
black hole region simply disappears as c — and the region (2) smoothly reaches the region 
(1). The region (3) can be approached by increasing c from region (2). The region where the 
black hole has a negative heat capacity (r^ < r/i < Tc) shrinks as c gets larger and only the 
most positive eigen mode remains in the region (3). The other two modes are getting close 
and become no longer normalizable modes when c reaches region (3) (c = z^^). We can see 
this tendency in FigJ7]in Appendix B. The limit to A = with c fixed is similarly discussed. 
As A approaches to zero, the value of positive eigenmodes approach to zero (oc —A) and these 
modes becomes non normalizable in the limit A ^ 0. This is clear in Figj8]in Appendix B. In 
order to discuss this transition, it is more proper to introduce a finite cavity. If we introduce 
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1 2 3 4 5 6 r/j 



Figure 4: The eigenvalue A in terms of Vh Figure 5: The eigenvalue A in terms of 
for c = 1 with A = -1. G [0, 1] for c = 0.01 with A = 0. = 

a finite cavity at r = for the case of zero cosmological constant, we have a large black hole 
Th tb with a same temperature and this large black hole has a positive heat capacity. Then 
positive eigenmodes A oc are expected and become non normalizable as oo. 

In the next section we discuss how we can understand these results. We are especially 
interested in the negative modes of small black holes with a positive heat capacity. It is 
interesting to notice that the number of negative modes in this region is two or one and 
there are two or one larger black holes with the same temperature for the negative or zero 
cosmological constant respectively. 

4 Local dynamics vs local thermodynamic vs global sta- 
bility 

In region (2), we find the discrepancy between the sign of heat capacity and the existence of 
negative modes. This leads us to reconsider the connection between the stability against small 
perturbations and the local and global thermodynamic stability of black hole. As mentioned 
in Introduction, the argument on the connection between the local thermodynamic stability 
and the stability against small perturbations is given in detail by Reall [3] , and they are shown 
to agree under some reasonable assumptions. 

It is good to start with looking at the counter examples [10] for the so-called correlated 
stability conjecture by Gubser and Mitra [5]. This conjecture is stated for a black string. 
Let us denote the direction to which a black string extends by z. The classical gravitational 
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S-mode perturbations having the wavefunction e^^^ along z satisfy the Lorentzian Einstein 
Lichnerowicz equation (13. ip of a black hole in a lower dimensions with a negative eigenvalue 
A = —k"^. Gubser and Mitra conjecture that a black string with translational symmetry and 
infinite extent shows a Gregory-Laflamme instability if and only if it has a local thermodynamic 
instability. In many cases, this conjecture is checked to be true. However the counter examples 
are known and those black strings show a Gregory-Lafiamme instability even they have the 
local thermodynamic stability. This discrepancy is understood as follow [10]. The black 
string is not uniquely identified by the conserved charges and has moduli along which the 
conserved charges do not change. Since local thermodynamic of black hole cannot discuss the 
stability along the moduli, the Gregory-Laflamme instability can exist along the moduli and 
this happens in the counter examples. 

Let us discuss whether a similar argument can apply for our numerical results where we 
found there exist two negative modes in a small black hole in region (2). First of all, there 
are two known black holes, e = 1 or —1 in (12.51) . However they have different values of 
the asymptotic AdS curvature in the infinity and they are not connected by a continuous 
parameter. A simple product of n dimensional Schwarzschild black hole in Gauss-Bonnet 
theory and 5*^ is not a solution of + 1 dimensional Gauss-Bonnet theory and then we 
cannot simply apply their argument on a black string instability here. If we can embed the 
Schwarzschild black hole metric in Gauss-Bonnet theory into a black string metric as a sub 
manifold, we may find a black string is not uniquely identified by the conserved charges. 

The non uniqueness of a black hole is discussed in [13] and the Reissner-Nordstrom black 
hole in four dimensions can be made classically unstable by developing a massive scalar hair 
under the mass and charge fixed. There is a massive scalar mode around the fiat space in 
Gauss-Bonnet theory \IA^, but the negative mode is a off-shell mode and we do not see a 
connection between them. 

Thus it is not clear that the existence of negative modes in a small black hole has a 
connection to the non uniqueness of a black hole. On the other hand, we have an observation 
that having two negative modes seems to have a relation with the fact that there are two larger 
black holes for a smallest black hole at a same temperature, which leads us to reconsider the 
global stability. (The number of negative mode is one for the case of zero cosmological constant 
and there is one larger black hole for a small black hole at a same temperature.) 

We then refer to the original interpretation of having a negative mode as the tunneling 
process which induces the decay of a metastable hot fiat space into nucleation of black holes 
in [TJ. At a temperature when the smallest black hole has two negative modes, the free energy 
of the smallest black hole is largest among that of three black holes and the free energy of the 
largest black hole is smallest. Then we may speculate that the large negative mode indicates 
the tunneling process of a hot AdS space into the nucleation of the largest black holes by 
passing through the smallest black hole solution. The small negative mode may indicate a 
more complicated tunneling process of a hot AdS space into the nucleation of the largest 
black holes by first passing through the small black hole solution and next passing through 
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Figure 6: The one loop effective potential around the AdS space (A), smallest (B) and middle 
size (C) black holes in region (2). The tunneling processes are shown by dotted lines. 



a middle size black hole (FigE])- We notice that we found positive eigenmodes around the 
middle size black hole. In the entropy comparison between a black string and a black hole 
in the context of Gregory-Laflamme instability, the instability mode is a classical (on-shell) 
mode and thus we have to compare the entropy with a fixed mass when we discuss about the 
global thermodynamic stability. However the negative modes are off-shell modes in quantum 
gravity and a tunneling process is discussed at a fixed temperature. 

In this section, we discussed some speculations from our numerical results, but everything 
is not conclusive at all. 



5 Summary and Discussions 

We studied non conformal negative modes around the Schwarzschild black hole in five dimen- 
sional Gauss-Bonnet theory with and without a negative cosmological constant. We found 
an unique negative mode when a black hole has a large radius and a negative heat capacity. 
The negative value goes to zero exactly when the sign of heat capacity changes its sign from 
negative to positive as the black hole radius becomes larger. This behavior is similar to the 
case of Einstein gravity and is consistent with the fact that the effect of Gauss-Bonnet term is 
negligible for a large black hole. On the other hand, we found negative modes still exist even 
when a small black hole has a positive heat capacity and the number is one or two for zero 
or non zero negative cosmological constant respectively. We do not fully understand why this 
happens and this result may indicate the black hole is not uniquely parametrized by conserved 
charges. Even we could refine the thermodynamic stability to be consistent with the existence 
of negative modes, it would be still difficult to predict the number of negative modes. Fol- 
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lowing the original argument in [T], we speculated that the negative modes may indicate the 
tunneling processes of the decay of a metastable hot AdS space into the nucleation of black 
holes. However these are speculations and a more further study is necessary. 

There is a report [I5] that a Schwarzschild black hole in five dimensional Gauss-Bonnet the- 
ory without the cosmological constant shows a classical instability against the perturbations 

which have arbitrary high angular momentum along for the range rh G [0, ^J'2(^^^\-~\/^)c] 
with c > 0. This instability is expected to exist, at least, for the small negative cosmological 
constant (— Ar^ <^ 1). This result is also counter intuitive since an angular momentum effec- 
tively plays as a mass and we do not expect there exist an instability by heavy massive modes. 
Therefore the black hole thermodjTiamics does not capture many gravitational dynamics. 

It is interesting to search a possibility to embed a Schwarzschild black hole in Gauss-Bonnet 
theory into a black string in a higher dimensions in order to discuss the correlated stability 
conjecture. In this direction it is interesting to study the Gregory-Lafiamme instability for a 
black string solution in Gauss-Bonnet theory constructed in [16j. 

The quantum effects become important when the spacetime curvature becomes strong and 
we expect a part of such quantum effects can be captured from studying the effects by higher 
dimensional operators like Gauss-Bonnet terms. It would be important to understand our 
results in terms of quantum gravity. 
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A Numerical analysis 

We explain the procedure of numerical analysis. We solve the equations of motion 



We can see most of the equations are solved with the above ansatz and finally have three 
nontrivial equations which is schematically written 



(A.l) 



with the ansatz 




(A.2) 




(A.3) 
(A.4) 
(A.5) 
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where Fi {i = 1,2,3) are functions of variables written inside the bracket and H[^{r) = 
dHtt{r) /dr etc. From counting the number of variables and equations, it is easy to see that 
there are two independent modes in these perturbations which we refer the 'traceless' and 
'trace' modes since they are traceless and trace modes in Einstein gravity limit. 

We solve the equations of motion from the horizon to the infinity. Because of the invariance 
under the rescaling (13.101) . we can take A = — 1 in the numerical analysis. The parameters are 
c, Th and A. We fix c and and change A to find eigenmodes by the shooting method. 

We analytically solve the equations of motion near the horizon up to e^, (e = r — r^), 
for Htt and !!„■ and up to for K{r) and compute the boundary condition at a small value 
e, i.e. Httijh + e), Hrrijh + e), Kir^ + e) and K'{rh + e). We took e = 0.001 for the cases 
A = -1, c = -0.1,0,0.01,1 and rn G [0,10] in FigHH We then solve the equations of 
motion numerically using a Fehlberg fourth-fifth order Runge-Kutta method with degree four 
interpolant in Maple 11 up to a large value rj in r, rj = 10^ in FigHMl 

We then fit the numerical values Hu{rf), Hrr{rf), K{rf) and K\rf) by the asymptotic 
solutions at infinity. We analytically compute the asymptotic solutions up to the second 
leading orders in r. There are four solutions since there are normalizable and non-normalizable 
solutions of 'traceless' and 'trace' modes. We fit the numerical values by using these all four 
solutions, since even we start the boundary condition near the horizon such that the 'trace' 
mode is set to be zero, there always exist numerical errors and the 'trace' modes are included 
in the numerical values. We them obtain the coefficient of the non-normalizable solution 
of 'traceless' mode in the numerical values and study when the coefficient becomes zero by 
changing A using the shooting method. 

We checked our numerical analysis has a good accuracy from various approaches. Our 
results approach to the one in Einstein gravity by taking c smaller values. We change the 
value of A and checked the invariance under the rescaling is satisfied with a good accuracy. 
We changed the value of e and rj and studied the results do not change. 

In the case of zero cosmological constant, the procedure is same except that rj = 100. 
Since the asymptotic solutions are exponential (oc e^^^~^^) for the 'traceless' mode and sine 
(oc sin(A/±2A/3)) for the 'trace' mode, the exponential growing is already manifest at = 
100. 

B Numerical results 

We show our numerical results. These values are used to plot the figures. We employ the 
shooting method and change A by 0.01 to find the eigenvalue. Thus for example the first value 
(0.8, —2.925) means the real value A is located between —2.92 and —2.93. 

(1) (c,A) - (-0.1,-1) : FigDD 
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(2) (c,A) = (0.01,-1) : Figl 
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(4) (c, A) = (0.01,0) : FigE 



rh 


0.0612 


0.062 


0.065 


0.07 


0.08 


0.1 


0.12 


0.14 


0.16 


X 


-0.015 


-0.065 


-0.315 


-0.835 


-2.095 


-4.595 


-6.415 


-7.465 


-7.925 




rh 


0.18 


0.2 


0.4 


0.6 


0.8 


1 


2 


4 


10 


X 


-7.965 


-7.765 


-3.835 


-1.975 


-1.175 


-0.765 


-0.195 


-0.055 


-0.0085 



(5) (c,A) = (0.1,-1) : FigH 
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(6) (c,A) = (0.01,-0.1) : Figl 
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Figure 7: A for (c,A) = (0.1,-1). The Figure 8: A for (c,A) = (0.01,-0.1). The 
dotted hue is A for (c. A) = (0.01, -1). dotted hue is A for (c. A) = (0.01, -1). 
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